The European Physical Journal E

Tumor spheroid mechanics arise from both the mechanical properties of individual cells and the adhesive interactions that organize them into tissues. The relative contribution of these two factors to the bulk mechanical behavior, however, remains difficult to disentangle experimentally. Here, we develop a computational model of micropipette aspiration to compare the mechanical response of isolated cells and multicellular spheroids within a common computational framework. By independently varying single-cell stiffness and cell-cell adhesion, we quantify their effects on aspiration dynamics, effective elastic modulus, and viscoelastic relaxation. Our results show that increasing single-cell stiffness substantially alters the mechanics of isolated cells but has limited influence on the effective elastic modulus of multicellular spheroids. In contrast, changes in cell-cell adhesion produce pronounced effects on spheroid effective elastic modulus. Nevertheless, both parameters increase the retardation time governing the transition from the initial elastic response to long-time viscous deformation. These findings suggest that multicellular elasticity is governed primarily by intercellular mechanical coupling, whereas the dynamical response to applied stress depends jointly on cell-scale mechanics and cell-cell adhesion.

Wound repair is an important mechanism to preserve tissue integrity in organisms after injury. However, why different tissues exhibit different mechanisms to repair wounds is a long-standing question that remains unanswered. In this work, we theoretically explore the role of the purse string, an actomyosin contractile cable used by tissues to close small wounds. Does the tissue 3D geometry influence the efficiency of the purse string in driving wound closure? Using a 3D biophysical model, we study in silico tissues with the same cell volumes but different aspect ratios, ranging from squamous to thick and tall tissues. The model predicts that taller cells are easily deformed by the purse string. In contrast, very squamous cells require a very strong purse string that might demand additional cellular mechanisms to close the gap. These findings establish a theoretical framework to predict the optimal biophysical mechanisms of wound healing in different tissues. Graphical abstractCells of different aspect ratios can be observed in a range of organisms with different function and mechanics. The wound healing efficiency of the purse string increases with the cell aspect ratio in our theoretical exploration. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/743165v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@d44ab0org.highwire.dtl.DTLVardef@1737cbaorg.highwire.dtl.DTLVardef@101b5d4org.highwire.dtl.DTLVardef@1487f26_HPS_FORMAT_FIGEXP M_FIG C_FIG

Electrostatic actuation is an emerging technology for generating tactile sensations on capacitive touchscreens through voltage-induced attractive forces between a fingertip and the surface. However, accurate control of electrostatic attraction during natural touchscreen interactions remains challenging because the applied normal force and sliding speed continuously vary, and their effects on the fingertip-screen contact and resulting actuation strength are not fully characterized. Here, we show how normal force and sliding speed systematically alter fingertip- screen contact area and electrical impedance, and use these measured changes to estimate electrostatic attraction during sliding. Contact area, interaction forces, and electrical impedance were measured simultaneously as participants slid their fingertips across an electrostatic surface under systematically varied normal forces and sliding speeds. These measurements revealed condition-dependent changes in fingertip contact, electrical interaction impedance, effective capacitance, derived effective gap thickness, and electrostatic attraction. We then incorporated these measured contact quantities into a physics-informed, data-driven model based on parallel-plate capacitor theory, in which effective capacitance, apparent contact area, and effective voltage determine the estimated electrostatic attraction. The resulting model links force- and speed-dependent changes in these quantities to electrostatic attraction while accounting for inter-participant variability through a participant-specific scaling factor. These findings provide experimentally grounded guidance for designing electrostatic surface-haptic feedback and future adaptive control strategies under realistic touch conditions.

Fingertip friction plays a central role in tactile exploration and object manipulation. During sliding, tangential force depends jointly on the real contact area and the interfacial shear stress, both of which can be influenced by sliding conditions. However, changes in fingertip friction are often interpreted primarily through changes in real contact area, whereas the accompanying changes in interfacial shear stress remain less well characterized. This gap is especially relevant for electrostatic surface haptic displays, which modulate fingertip friction by applying a voltage between the finger and the touch surface. Here, we experimentally quantify the mean interfacial shear stress of a sliding fingertip on an electrostatically actuated touchscreen using simultaneous measurements of tangential force and optically resolved real contact area. Ten participants performed sliding trials across three speeds and three normal forces with and without electrostatic actuation. Interfacial shear stress increased with speed and decreased with normal force; in both cases, these trends arose because real contact area varied more strongly than tangential force. Electrostatic actuation further reduced interfacial shear stress, as increasing voltage produced a larger increase in real contact area than in tangential force. These findings show that interfacial shear stress varies systematically with sliding conditions and electrostatic actuation, clarifying how changes in real contact area and interfacial shear stress combine to shape fingertip-surface friction.

Epithelial tissues continuously remodel their curvature during morphogenesis, homeostasis, regeneration, and disease, yet experimental access to time-varying curvature remains limited. Here, we introduce CurvoChip, a pneumatically actuated microsystem that reversibly deforms confluent epithelial monolayers cultured on a 20-m elastic membrane into concave or convex geometries. The device operates either in a standard incubator or on a microscope stage and provides programmable control over pressure amplitude, direction, and cycling. Analytical scaling, finite-element simulations, and confocal profilometry establish predictable membrane deformation across the operating range, whereas cycling between -400 and +400 mbar for 120 cycles produces stable deflection without detectable drift or residual deformation. We further implement a three-dimensional surface-reconstruction and segmentation workflow to quantify cell and nuclear morphology on curved monolayers. Acute curvature induction produces a marked polarity-dependent response: convex deformation causes greater cell spreading and epithelial thinning than concave deformation, while nuclear projected area, thickness, and volume change in a direction- and position-dependent manner. These results show that epithelial architecture is sensitive not only to curvature magnitude but also to its orientation relative to the apico-basal axis. CurvoChip therefore provides an accessible platform for dissecting how epithelial tissues integrate dynamic geometric cues.

While geometric constraints shape tissue development, quantifying the resulting growth dynamics remains a central challenge in tissue engineering. Conventional methods often struggle to capture multi-scale kinetics without complex labeling or difficult single-cell tracking. Here, we analyze geometrically controlled growth of microtissues derived from human dermal fibroblasts using time-resolved, label-free brightfield microscopy, combined with optical flow and semi-automated deep learning mitosis detection. By extracting multi-scale flow fields and integrating them with tissue segmentation, we quantify directional tissue dynamics, separating flow into parallel and normal components relative to the local tissue contour. Applying this framework, we contrast the quiescent tissue interior with the advancing growth front where localized dynamics and cell proliferation drive expansion. Our results demonstrate that, compared to the bulk, the growth front exhibits higher fluctuations parallel to the tissue contour, positive mean normal flow, and significantly increased mitotic activity. Furthermore, evaluating flow divergence around mitotic events reveals distinct spatial behaviors: with the onset of mitosis, a contraction and subsequent expansion occurs in the vicinity of the dividing cells. Beyond the immediate cellular neighborhood, the broader regional dynamics remain consistent before and after mitosis onset, with net tissue expansion in proximity to the growth front and contraction within the tissue interior. By extracting continuous kinetic data from easily accessible, label-free brightfield imaging, this approach serves as a non-invasive, complementary tool for evaluating in vitro tissue morphogenesis and growth dynamics. This analytical framework can be expanded to study locally resolved tissue morphogenesis and growth kinetics in other microsystems, ranging from embryos to organoids. Statement of SignificanceUnderstanding how localized cellular forces drive tissue growth is critical for mechanobiology. However, mapping these dynamics traditionally requires complex, invasive fluorescent labeling. We present an accessible, label-free computational framework combining optical flow and deep learning-based mitosis detection to quantify continuous tissue kinematics directly from standard brightfield microscopy. Applying this to 3D microtissues, we reveal a distinct spatial coupling between cell division, local mechanical fluctuations, and directed tissue expansion at the active growth front. This non-invasive approach bridges the gap between single-cell mechanics and macroscopic morphogenesis, offering a versatile tool to monitor complex in vitro model systems-like organoids and bioengineered tissues-without disrupting their native state.

Mechanical properties of biological tissues, driven by passive and active forces, play a vital role in several processes ranging from development to cancer metastasis. However, the dynamical responses of cells in tissues, subject to mechanical deformations such as shear and the associated rheological properties, are not well characterized. Here, we use three-dimensional agent-based models for normal and cancer tissues to investigate their responses to simple shear as a function of cell stiffness and stochastic active forces. In the normal epithelium, with uniform strength of active force, the yield stress as a function of shear rate follows the Herschel-Bulkley form over a range of cell volume fraction. Strikingly, the shear rate dependence and the elasticity-dependent changes in the yield stress fall on master curves upon suitable scaling. To model cancer-like behavior, a certain fraction (Np) of cells was chosen to have enhanced activity and decreased stiffness. As Np increases, the extent of collective cell movement decreases, transitioning from affine (collective) to non-affine (individualistic) movement, a finding that is in accord with imaging experiments. Simulations of a model of a stiff solid tumor, with radius Rs embedded in normal tissue, show that as Rs increases, the yield stress increases. Interestingly, the cells migrate collectively as Rs increases. A Gaussian Mixture Model (GMM) and a mean field theory quantitatively account for the simulation as well as experimental results on cancerous, non-cancerous, and a mixture of these two types. The combined theoretical and experimental study establishes that heterogeneity in stiffness and activity determines non-affine movements in normal and cancer tissues.

Polyglutamine (polyQ) aggregation plays a central role in several neurodegenerative diseases, including Huntington's disease. DNAJB6b, a molecular chaperone involved in protein quality control, is known to efficiently suppress polyQ aggregation, but its anti-aggregation mechanism remains unclear. In this work we investigate the interaction between DNAJB6b and the polyQ region (Q48) of mutant Huntingtin Exon 1 (mHttEx1) using a custom-built coarse-grained molecular dynamics model. The model incorporates a two-bead-per-amino-acid representation with hydrogen bonding (termed 2BPA-HB), and is calibrated against all-atom molecular dynamics data in terms of geometry, hydrophobicity, and hydrogen bonding. The model reproduces the tertiary structure of DNAJB6b and its interactions with Q48, and reveals an inverse correlation between DNAJB6b concentration and Q48 aggregation propensity. Our simulations show that DNAJB6b co-condensates with polyQ molecules, thereby shielding the polyQ from forming the intermolecular hydrogen bonds necessary for amyloid formation. The 2BPA-HB CGMD model en- ables efficient exploration of DNAJB6b conformations, supporting future studies of chaperone-mediated aggregation suppression and therapeutic development.

Cardiac blood flow is a regulator of several important developmental and remodelling processes in the heart, including through fluid shear forces sensed by the endothelial cells lining the heart. However, optically mapping these flow fields in the complex 3D geometry of the heart is challenging even in transparent animal models such as the zebrafish. One of the main challenges is the difficulty in measuring the out-of-plane (axial) velocity component, preventing accurate mapping of the complete 3-component-3-dimension (3C-3D) blood flow velocity field; image-based techniques such as microscopic particle image velocimetry ({micro}PIV) traditionally only provide the in-plane flow components. Here we present a computational approach to achieve full time-varying 3C-3D blood flow vector mapping using a standard selective plane illumination microscope (SPIM), based on robust cardiac phase assignment, precise measurement-driven registration of sequentially acquired z-stacks, and PIV data fusion from multiple sample orientations. Our approach holds the key to understanding the complex dynamic flow fields within the developing heart, and their role in shaping cardiac development.

A classical study found no excitation transfer when isolated cardiomyocytes were placed side by side, whereas a recent paper reported action potential transfer in carrdiomyocytes placed end to end. We use nanoscale numerical simulations based on the full Poisson-Nernst-Planck equations to investigate whether these apparently opposing observations can be explained by the different geometrical configurations. The computations show that in the end-to-end configuration, ephaptic coupling occurs when the intercellular cleft is sufficiently narrow and a sufficiently large fraction of the sodium channels is localized at the intercalated disc. Coupling is strengthened when the sodium channels are concentrated in fewer clusters and when ionic diffusion within the cleft is reduced. Under these conditions, excitation transfer occurs on a timescale consistent with rapid cell-to-cell activation. Conduction depends biphasically on cleft width and terminates abruptly beyond a critical width. Localization of potassium channels at the intercalated disc has only a moderate effect, whereas gap junctions substantially improve conduction and reduce the relative contribution of ephaptic coupling. In the side-by-side configuration, excitation transfer does not occur under physiological conditions and requires highly flattened cells, minimal separation, and unrealistically strong sodium-channel clustering. The different outcomes of the side-by-side and end-to-end experiments can therefore be explained by the fundamentally different geometrical conditions for ephaptic coupling.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

Classical neuronal cable theory relies on quasi-static electric field approximations and neglects magnetic induction, Lorentz force coupling, and transient electromagnetic currents, limiting its ability to fully characterize action potential propagation within geometrically branched axons and dendrites. This work develops a coupled Maxwell-electromagnetic cable framework by integrating finite-difference time-domain (FDTD) solutions of Maxwells equations with extended Hodgkin-Huxley and Fitzhugh-Nagumo membrane dynamics, incorporating magnetic gating perturbations, electromagnetic trans-membrane currents IEM, and nanoscale quantum corrections for thin neural segments. Controlled propagation experiments are designed to quantify deviations from standard cable predictions across asymmetric and symmetric axonal bifurcation geometries. Numerical results demonstrate that inductive magnetic effects lower the critical branch radius for junction conduction failure and break symmetric action potential invasion in geometrically identical child branches under external transverse magnetic fields. An electromagnetic corrected geometric ratio GREM is proposed to revise impedance-matching conditions at branch points, accounting for size-dependent axial current imbalance induced by magnetic and displacement currents. Parent axon conduction velocity deviates substantially from the canonical [Formula] scaling law when electromagnetic feedback and quantum charge distributions are included, triggering early signal blockage at large cable diameters. Collectively, this study establishes that quasi-static cable models underestimate electromagnetic corrections to propagation speed, waveform shape, and bifurcation transmission fidelity; the coupled Maxwell-cable framework provides a comprehensive multi-physics tool for modeling electrodynamic signal behavior in complex neuronal architectures.

Quantification of microbial growth inhibition is central to assays ranging from antibiotic susceptibility of bacteria to sensitivity of yeasts to antifungal therapeutics. Classical analysis approaches derive from zone-of-inhibition (termed halo) formats using filter paper discs, spanning methods from laser detection to machine learning. However, these tools struggle with non-uniform halos, fail to account for lawn density variability despite its experimental influence, and lack accessible, reproducible code. Here, we present Halo Unbiased Measurement of growth Inhibition (HaloUMI); an open-source Python graphical user interface for automated, high-throughput analysis of lawn-based microbial assays. HaloUMI integrates robust image processing with physics-informed models to quantify inhibition zones irrespective of shape, enabling accurate segmentation of uniform and irregular halo phenotypes. This analysis pipeline incorporates the critical correction for spatial heterogeneity in lawn density, improving reproducibility across experimental conditions. The software enhances usability without sacrificing precision, allowing rapid batch processing and intuitive parameter control. HaloUMI can be applied to multiple assay types, including yeast toxin halo, microbial mating, and conventional filter paper disc assays. It yields high-precision measurement of halo size and morphology, with improved consistency compared to standard thresholding and circular fitting. By combining accessibility, flexibility, and biophysical modelling, HaloUMI provides a quantitative framework for irregularly shaped halos of lawns of varying growth potential, enabling generalisable analysis of broad microbial interactions. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=84 SRC="FIGDIR/small/743694v1_ufig1.gif" ALT="Figure 1"> View larger version (20K): org.highwire.dtl.DTLVardef@1bd7c88org.highwire.dtl.DTLVardef@13ac9b6org.highwire.dtl.DTLVardef@9108e1org.highwire.dtl.DTLVardef@1de06bd_HPS_FORMAT_FIGEXP M_FIG C_FIG

Based on our experimental observation of activation state oscillations of different frequencies used to communicate information by the master human stress response regulator protein p38 MAPK 1, we develop a simple graphical approach for understanding and predicting the behavior of complex biological networks acting upon signals carrying information as different frequency waves of chemicals. This approach uses the same techniques used for analyzing and understanding information transmission using waves of electrical currents and fields used in electrical alternating current (AC) circuits. We show how biological components can be organized to behave as standard components found in electronic telecommunications circuits. Finally, we demonstrate how such components can be organized into complex biological signaling cascades whose behavior can be qualitatively and quantitatively understood, and whose output resembles that experimentally observed in p38.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

The study of tissue dynamics has been stimulated during the last decades thanks to the use of quantitative descriptions, with the development of several theoretical and computational frameworks, many of them revolving around the notion of reaction-diffusion systems, eventually with additional structure variables beyond time and space. The use of structure variables can accommodate phenotypic traits. In this work, we study a family of competition models, where a given population depends on a resource (e.g. oxygen) and several populations are competing for it. Our quantitative description incorporates phenotypic traits and heterogeneity at the level of cell cycle variations, which influence replication rates via oxygen consumption. This enables us to replicate the fitness of specific subpopulations to environmental conditions (e.g. oxygen shortage or external influences). Using numerical simulations, we show that such models display dynamical pattern formation in the form of coupled travelling wave profiles that expand or retreat at the same wave speed. The full theoretical analysis of such dynamics is quite involved; to circumvent this difficulty, we introduce a quasi-stationary approximation for the resource dynamics. We find that this approximation can reproduce the overall behaviour very accurately, with the additional benefit of allowing theoretical treatment of the reduced model. In this way, we provide estimates on the wave speed which are numerically shown to be robust across a wide range of macroscopic parameters of the full model. The wave speeds are thus found to depend strongly on the proliferation rate of the fittest population, resembling a winner-takes-all dynamics.

Recent experiments show that protein condensates sitting on opposite surfaces of a flat lipid membrane move together and prefer to overlap, even though they cannot touch each other. This points to an indirect, membrane-mediated interaction. Two mechanisms could be responsible: a curvature-induced interaction, which is energetic in origin, and a fluctuation-induced interaction, which is entropic. Here we study both with coarse-grained molecular dynamics simulations, using Cookes implicit-solvent lipid model together with a generic bead-spring polymer model for the condensate. We compute the potential of mean force between two condensates across the membrane. For condensates of the same size, full overlap is unfavorable, and the pair instead settles into a partially overlapping state that bends the membrane into an S-like shape. When the two condensates differ strongly in size, full overlap becomes favorable. We explain this with a simple geometric picture. The condensate wets the membrane as a thin film and imposes curvature only along its rim, while membrane tension flattens the membrane under its interior. The resulting ring of curvature can trap a smaller condensate on the opposite side. We also compare the bending undulations and the effective bending modulus of a bare membrane, a membrane with one condensate, and a membrane with condensates on both sides. A wetting condensate suppresses the undulation modes and stiffens the membrane, but whether this makes overlap entropically favorable remains inconclusive. Our results indicate that the coupling is driven mainly by curvature, and that it depends on the wetting mechanism and on the membrane tension.

In the active loop extrusion model, the cohesin protein complex creates chromatin loops in eukaryotic cells. Extrusion maintains topologically associated domains (TADs), which are contiguous segments of chromatin that preferentially colocalize in space and are typically bounded by CTCF proteins that pause cohesin translocation. Here, we model active loop extrusion with hybrid molecular dynamics - Monte Carlo simulations in entangled flexible linear polymer melts. Intra-chain contact probabilities of polymers with active loop extrusion are enhanced compared to their equilibrium, passive counterparts. Extrusion causes the size of chain segments to be much smaller than in passive melts. While the overlap parameter in passive melts without extrusion monotonically increases with segment length, it is nonmonotonic in active melts and on the order of unity within the parameters of this study. Active loop extrusion suppresses contacts between TADs in favor of intra-TAD contacts. Reduction of overlaps between chain segments dilutes entanglements in active melts. Depending on parameters, active extrusion without TADs may induce more compact conformations than with TADs, due in part to fractal loopy globule-like dynamics. This work suggests that active loop extrusion reduces overlaps between TADs, contributing to effective gene regulation by cis-regulatory elements.

This study aimed to (1) develop an elasticity framework for the sprint force-velocity (F-V) relationship and (2) examine how maximal force (F_{0}), maximal velocity (v_{0}), and sprint distance modulate the four derived elasticity metrics, and (3) explore these elasticity metrics' interrelation. After modelling the F-V relationship differential equation, four elasticity metrics were defined as force elasticity (F_{e}), the elasticity of sprint time to F_{0}; velocity elasticity (v_{e}), the elasticity of sprint time to v_{0}; the force-velocity elasticity norm {(\mathrm{F}-\mathrm{V}}_{\mathrm{EN}}=\sqrt{F_{e}^{2}+v_{e}^{2}}), capturing the combined sprint time sensitivity to proportional changes in F_{0} and v_{0}; and the force-velocity elasticity ratio {(\mathrm{F}-\mathrm{V}}_{\mathrm{ER}}=F_{e}{\div v}_{e}), indicating which variable dominates the sprint time response. Model simulations showed that F_{e} decreased with rising F_{0} and increased with rising v_{0}, while v_{e} showed the opposite pattern. With increasing sprint distance, F_{e} decreased and v_{e} increased. Given its negligible effect on sprint time, ignoring air resistance yields a conservation law (2F_{e}+v_{e}\equiv 1), indicating that a gain in one elasticity metric necessarily diminishes the other in a fixed proportion. This framework also identifies a valley distance (d_{valley}) at {\mathrm{F}-\mathrm{V}}_{\mathrm{ER}}=2, where {\mathrm{F}-\mathrm{V}}_{\mathrm{EN}} is minimized (\sqrt{0.2}) and sprint time is least responsive to changes in F-V relationship variables. Empirical data confirmed that the two theoretical laws still hold approximately when air resistance is considered. By linking changes in F_{0} and v_{0} to sprint time across different distances, the elasticity framework provides a quantitative basis for estimating the theoretical sprint time response to documented changes in F-V relationship variables.

Synapsin-1 is a multidomain neuronal protein containing extensive intrinsically disordered regions and is a key component of synaptic-vesicle condensates. Direct residue-level simulation of the collective organization of thousands of synapsin molecules remains computationally demanding. Here, we develop a coarse-grained description that connects residue-level CALVADOS 3 simulations to a one-particle-per-protein model. A potential of mean force between two synapsin molecules is obtained by umbrella sampling and represented by an isotropic effective interaction containing a short-range attractive region and a weak outer repulsive contribution. We compare two treatments of this interaction that differ only in the retention of the outer tail. Langevin dynamics simulations of effective proteins show aggregation upon cooling and compression in both models, but with markedly different collective organization. The shorter-ranged model progressively coarsens toward a single dense domain, whereas retaining the outer repulsive contribution favors the persistence of multiple mesoscale aggregates. The two models also display distinct relationships between aggregate size and particle mobility at low temperature. These results show that weak features of an effective protein-protein interaction can have pronounced consequences for collective synapsin organization at mesoscopic scales. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=95 SRC="FIGDIR/small/745062v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@854288org.highwire.dtl.DTLVardef@d2fc52org.highwire.dtl.DTLVardef@1b37d5aorg.highwire.dtl.DTLVardef@eac583_HPS_FORMAT_FIGEXP M_FIG C_FIG